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Abstract. Cosmologies with running cosmological term p\ and gravitational Newton's 
coupling G may naturally be expected if the evolution of the universe can ultimately 
be derived from the first principles of Quantum Field Theory or String Theory. For 
example, if matter is conserved and the vacuum energy density varies quadratically 
with the expansion rate as p\{H) = no -|- n2H^, with no / (a possibility that has 
been advocated in the literature within the QFT framework), it can be shown that 
G must vary logarithmically (hence very slowly) with H. In this paper, we derive the 
general cosmological perturbation equations for models with variable G and p\ in which 
the fluctuations 6G and 6p\ are explicitly included. We demonstrate that, if matter is 
covariantly conserved, the late growth of matter density perturbations is independent 
of the wavenumber k. Furthermore, if p\ is negligible at high redshifts and G varies 
slowly, we find that these cosmologies produce a matter power spectrum with the same 
shape as that of the ACDM model, thus predicting the same basic features on structure 
formation. Despite this shape indistinguishability, the free parameters of the variable 
G and p\ models can still be effectively constrained from the observational bounds on 
the spectrum amplitude. 

PACS: 95.36.-hx 04.62. -hv ll.lO.Hi 



1 Introduction 



The current "standard model" (or "concordance model" ) of our universe, being an homogeneous 
and isotropic FLRW cosmological model, consists of a remarkably small number of ingredients, to 
wit: matter, radiation and a cosmological constant (CC) term, A. The first two ingredients are 
dynamical and vary (decrease fast) with the cosmic time t whereas the third. A, remains strictly 
constant. After many years of theoretical insight (cf. e.g. the reviews [Tl[2] and references therein) , 
the situation of the original FLRW cosmological models has not changed much, in the sense that 
we have not been able to make any fundamental advance in the comprehension of the relationship 
between the matter energy density and the CC. Still, we have performed a major accomplishment 
at the phenomenological level by simultaneously fitting the modern independent data sets emerging 
from LSS galaxy surveys, supernova luminosities and the cosmic microwave background (CMB) 
anisotropics [3HZ]- On the basis of this successful fit, in which A enters as a free parameter, one 
claims that a non-vanishing and positive cosmological constant has been measured. Nevertheless, 
we still don't know what is the true meaning of the fitted parameter A. Assuming that Newton's 
gravitational coupling G is strictly constant, the combined set of observational data determine the 
value 

PA = ^=^(2.3xlO-3ey)^ (1) 

which we interpret as the vacuum energy density. It corresponds to 17^ = pa/Pc — 0-7 when the 
PA density is normalized with respect to the current critical density p'^ = {sVh x 10~^ eV)^ - for 
a reduced Hubble constant value of /i ~ 0.70. Assuming (in the light of the same observational 
data) that the universe is spatially flat, this means that the current matter density normalized to 
the critical density is = p^/Pc — 

All the efforts to deduce the value of the energy density ([I]) - a very small one for all particle 
physics standards, except if a very light neutrino mass is the only particle involved [8] - have 
failed up to now. The main stumbling block to a solution is the fact that any approach based on 
the fundamental principles of quantum field theory (QFT) or string theory lead to some explicit 
or implicit form of severe fine-tuning among the parameters of the theory. The reason for this is 
that these theoretical descriptions are plagued by large hierarchies of energy scales associated to 
the existence of many possible vacua. 

This difficulty became clear already from the first attempts to treat the dark energy component 
as a dynamical scalar field [9]. The idea was to let such a field select automatically the vacuum 
state in a dynamical way, especially one with zero value of the energy density. More recently, this 
approach took the popular form of a "quintessence" field slowly rolling down its potential and has 
adopted many different faces [10] - for a review, see [2] . But the situation at present is even harder 
because the quintessence field ~ or, more generally, a dark energy (DE) field - should be able to 
naturally choose not zero but the very small number ([1]) as its ground state. 

Despite the fact that a working dynamical mechanism able to choose the correct vacuum state 
has to be found yet, another important motivation for the quintessence models is that they aim at 
explaining the puzzling coincidence between the present value of pA and the value of the matter 
density, p^. In other words, why Q^/Q^^ = 0(1)? Arguably, a dynamical DE should be a starting 
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point to understand this puzzle. Detailed analyses of these models exist in the literature, including 
their confrontation with the data |11H14| . 

On the other hand, the possibility that the cosmological term is a running quantity, which could 
be sensitive to the quantum matter effects, seems a more appealing ansatz, as it could provide 
an interface between QFT and cosmology [S1[T2]. This fundamental possibility has been recently 
emphasized in [16]- for a review, see [17]. Actually, the analysis of the various observational data 
show (see Ref. [18l[19]) that a wide class of dynamical CC models models are indeed able to fit 
the combined observations to a level of accuracy comparable to the standard ACDM model. In 
some cases, the dynamical nature of A allows these models to provide a clue to the coincidence 
problem |20t l21j. and maybe eventually to the full cosmological constant problem [22) . 

In general, in this kind of dynamical CC scenarios we have A = A(,^) where = ^(t) is 
a cosmological variable that evolves with time, and therefore ultimately A = A{t) is also time 
varying. Originally, these models were purely phenomenological, with no relation to QFT [23j . 
Nowadays we know that not all models of this kind are allowed, and the fact that the observational 
data can discriminate which of them are good candidates and which are not so good gives some 
insight into the function ^ = ^(t) |18] . A particularly interesting class of variable CC models 
are those in which the gravitational coupling G changes very slowly (logarithmically) with the 
expansion of the universe [24p25j. This scenario is possible e.g. if matter is covariantly conserved. 

In practice, ^ could be the expansion rate H = H{t), the scale factor a = a{t), the matter 
energy density, etc. In this paper, we shall assume that ^ = ff, similarly as it was done for models 
with a variable A interacting with matter |15 y 26 [ [27]. With this ansatz, we shall study the impact 
on the structure formation, showing that these models predict the same shape for the matter 
power spectrum as the ACDM model, a fact that would not apply e.g. for variable CC models in 
which the CC decays into matter. We find this feature remarkable and we shall explore its possible 
consequences and also some possible phenomenological tests of this kind of models. 

The paper is organized as follows. In the next section, we classify the possible scenarios with 
variable p\ and G. In section 3, we present the general set of coupled perturbation equations 
involving 5p\ and 6G, showing that the late growth of matter fluctuations becomes independent 
of the scale k. A class of models with variable cosmological term as a series function of the Hubble 
rate is introduced in section 4. In section 5, we concentrate on a particular model in this class, 
which is well motivated from the QFT point of view, and analyze the constraints imposed by 
primordial nucleosynthesis and structure formation. In the last section, we present and discuss 
our conclusions. Finally, an appendix is included at the end where we discuss gauge issues. 

2 Generic models with variable cosmological parameters 

We start from Einstein's equations in the presence of the cosmological constant term, 

Rfiu - ^gfiuR = 8irG {T^y + g^^ p^) , (2) 

where T^y is the ordinary energy-momentum tensor associated to isotropic matter and radiation, 
and p\ represents the vacuum energy density associated to the CC. Let us now contemplate the 
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possibility that G = G{t) and A = A(t) can be both functions of the cosmic time within the context 
of the FLRW cosmology. It should be clear that the very precise measurements of G existing in the 
literature refer only to distances within the solar system and astrophysical systems. In cosmology 
these scales are immersed into much larger scales (galaxies and clusters of galaxies) which are 
treated as point-like (and referred to as "fundamental observers", comoving with the cosmic fluid). 
Therefore, the variations of G at the cosmological level could only be observed at much larger 
distances where at the moment we have never had the possibility to make direct experiments. In 
practice, the potential variation of G = G{t) and A = A{t) should be expressed in terms of a 
possible cosmological redshift dependence of these functions, G = G{z) and A = A(z). 

Let us consider the various possible scenarios for variable cosmological parameters that appear 
when we solve Einstein's equations ^ in the flat FLRW metric, ds^ = dt^ — a'^{t)dx.'^. To start 
with, one finds Friedmann's equation with non- vanishing p\, which provides Hubble's expansion 
rate H = a/a (a = da/dt) as a function of the matter and vacuum energy densities: 

H^=^-^{p,n+PK). (3) 
On the other hand, the general Bianchi identity of the Einstein tensor in ([2]) leads to 

[G{T^u+9^.uPK)]=Q. (4) 

Using the FLRW metric explicitly, the last equation results into the following "mixed" local con- 
servation law: 

^ [G(p„ + pa)] + 3GF(p„+p„) = 0. (5) 

If /9A 7^ 0, then pm is not generally conserved as there may be transfer of energy from matter- 
radiation into the variable pA or vice versa (including a possible contribution from a variable G, 
if G ^ 0). Thus this law mixes, in general, the matter-radiation energy density with the vacuum 
energy p\. However, the following particular scenarios are possible: 

• \) G =const. and /ja =const. This is the standard case of ACDM cosmology, implying the 
local covariant conservation law of matter-radiation: 

Pm + '^H {pm+Pm) (6) 

• ii) G =const and p\ 7^ 0, in which case Eq.jS]) leads to the mixed conservation law 

Pk + Pm + ^H {pm + Pm) = 0; (7) 

• iii) G / and pA =const, implying G(pm + pa) + G\pm + 'iH^Pm + Pm)] = 0; 

• iv) G 7^ and p^ 7^ 0. In this case, if we assume the standard local covariant conservation 
of matter-radiation, i.e Eq. ([6]), it is easy to see that Eq. ([5]) boils down to 

(Prn, + PA)G' + GpA = 0. (8) 



4 



Notice that only in cases i) and iv) matter is covariantly self-conserved, meaning that matter 
evolves according to the solution of Eq. ([6]): 

p^(a) = p^a-°- =p^(l + z)"'", a^ = 3(l + w™), (9) 

where is the current matter density and ujm = Pm/ Pm = 0, 1/3 are the equation of state (EOS) 
parameters for cold and relativistic matter, respectively. We have expressed the result Q in terms 
of the scale factor a = a{t) and the cosmological redshift z = (1 — o)/a. 

In cases ii) and iii), instead, matter is not conserved (if one of the two parameters p^ or G 
indeed is to be variable, respectively). Explicit cosmological models with variable parameters as 
in case ii) have been investigated in detail in |15 y i9 1 [27]. Cosmic perturbations of this model have 
been considered in |28fl32j . Case iii) has been studied at the background level in [33]. Finally, the 
background evolution of case iv) been studied in d.ile.ent c^^te.ts ■„ Ret. In the 

next section, we shall address the calculation of cosmic perturbations in a general model where 
p\ and G can evolve with the expansion, and we will then specialize the set of equations for the 
type iv) model in which matter is covariantly conserved. Only in section [J] we will further narrow 
down the obtained set of perturbation equations for a concrete model with running cosmological 
parameters |24j . 



3 Perturbations with variable A and G 

For the analysis of the cosmic perturbations in the general running pA and G model iv) of 
the previous section, we have to perturb all parts of Einstein's equations that may evolve with 
time; namely the metric, the energy-momentum tensor for both matter and vacuum, and finally 
we must perturb also the gravitational constant. Einstein's equations ([2]) can be conveniently cast 
as follows: 

R^.. = SttG (t^, - ^g^uT^^ . (10) 

As a background metric, we use the flat FLRW metric: 

ds"^ = g^^dxf'dx'' = dt^ - a^{t) 6ijdx'dx^ . (11) 

In perturbing it, Qfj^^ — >• Qf^i, + h^^, we adopt here the synchronous gauged (^oo = ^Oi = 0). The 
total energy momentum tensor of the cosmological fluid can be written as the sum of the matter 
and vacuum contributions: 

T.' = (matter) + = -PT^i + {PT+PtWU. , (12) 

where 

PT = Pm + PA 

PT = Pm+Pk = UJmPm " PA , (13) 



^Some potential astrophysical implications of scenario iv) have been addressed first in [24] and recently in [34] 
^See the Appendix for an alternative calculation in the Newtonian gauge. 
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(notice that wa = Pa/ PA = — 1 even if the CC is running). As for the density and pressure, we 
introduce the perturbations in the usual form, except that we have to include also the perturbation 
for the vacuum energy density because, in the present context, A is an evolving variable. Thus, 
we have pj — )• pj + 6pi , pi — )• pi + Spi, where i = m,A. As warned above, if we allow for a variable 
gravitational coupling G, we must also consider perturbations for it: 

G^G + SG. (14) 

Obviously, the perturbations 6G and 5pA are the two most distinguished dynamical components 
of the present cosmic perturbation analysis, as they are both absent in the ACDM model. 

Finally, we have to perturb the 4- velocity of the matter particles, — t- + 5U^. For an 
observer moving with the fluid (i.e. a comoving observer) it reads: = {1,5W). 

We are now ready to find the perturbed equations of motion for this running model. As 
the fundamental equations describing the perturbations for our model, we will take: i) the 00 
component of the Einstein equations (jlOp : and ii) the generalized Bianchi identity (jl]), which splits 
into energy conservation (notice that G 7^ in the present framework) 

(GTo^) = {GT,n + G [r^^To- - F^qT^^ = , (15) 

and momentum conservation: 

{GTt) = d, (GTt) + G [T^^^rr - T';,,T^] = . (16) 

As for the 00 component of the Einstein equations, and taking into account that G gets perturbed 
as in (I14p . we have: 

0th order: — 3— = AttG {px + ^Pt) , 
a 

(17) 

1st order: h + 2Hh = Svr [G {6pT + 36pt) + dG {pr + ^Pt)\ , 
d ( /ifcfc \ 

where we have defined h = — \ — ^ ) and repeated Latin indices are understood to be summed over 



the values 1, 2, 3. As energy-density components, we have matter and the running cosmological 
constant, with EOS parameters indicated in (jl3p . and thus the perturbed equation takes on the 
form: 

h + 2Hh = 87r{G[(l + 3a;„)5p„ + 6pA + 3Spa] + SG[{1 + 3uj,n)Pm - 2pa] } . (18) 

Let us now work out the equation of local covariant conservation of the energy, Eq. (jl5p . After a 
straightforward calculation, the final equations read as follows. On the one hand, 

0th order: G{p^, + pa) + G{pm + Pa) + 3GHpm{l + Um) = ■ (19) 

If we assume local covariant conservation of matter, ^ fi^Q (^^g^ncr) ~ ^' ^'^^ ®' ^® 
(fT9|) indeed reduces to Eq. ([8]), as it should be. On the other hand. 



1st order: G 



5pm + SpA + pm(l + COm) [^^ ~ ^ 3ij[(l + U}m)Spm + {SpA + Spa)\ + 

+G{6pm + Spa) + [pm + PA + 3H{1 + uj^)pm] 6G + {p^ + pa)SG = . (20) 



6 



In the previous equation, we have introduced the variable 



(21) 



which represents the perturbation on the matter velocity gradient, and we have used 5U^ = 0. If 
we invoke once more the local covariant conservation of matter, both at 0th order - see Eq. ([6]) - 
and at 1st order, 

h 



Spm + /Om(l + Wm) 



3i?(l +Wm)5pm = 0, 



(22) 



equation (|20p can be further reduced to the simpler form: 

G [SpK + m{5pK + 6pk)] + G{5pr,^ + 6pa) + paSG + (p„ + pa)5G = . (23) 



Finally, let us work out the equation of local covariant conservation of momentum, Eq. (jl6p . In this 
instance, there are no 0th order terms (in the background, the momentum conservation is trivial). 
The perturbed result reads 



a^dt[Gp{l + uj)6W] + 5ahGp{l + uj)6U' + Gdi{5p) + pdi{6G) = , 



(24) 



where it is understood that we should sum over all the energy components, in this case, matter 
and A. By applying the operator di and transforming to the Fourier space, we get: 

fc2 



GpO + G{pe + pe + 5Hp9) 



[G6p + uop6G] = , 



(25) 



which, after summing over matter and CC, becomes: 

r 1 k'^ 

(1 + UJm) GpmO + GipmO + PmO + 5H Pr, 



[G{6pA + U)mSpm) + {(^mPm " Pa)SG] . (26) 



The part corresponding to matter conservation (^V^T^^^^^^^^^^ = Oj is 

r • 1 /c^ 

(1 + OJrn) PmO + PmO + 5H PmO = -^UJ^Spm , 



SO that Eq. (|26p simplifies as follows: 



(1 + UJm)GpmO = [G5pA + {l^mPm " Pk)5G] 



(27) 



(28) 



Summarizing, our final set of equations is given by p8|) . (p2]) . (p3|) . (p7|) and ([28]) . It is particularly 
relevant for our purposes to rewrite these equations in the matter dominated epoch (w^ = 0) and 
assuming adiabatic perturbations for the CC {5pA = —^Pa)- In a nutshell, we find: 



h + 2Hh = Svr [pr,^ - 2pa] 5G + BvrG [6pm - 26 pa] 



Spm + Pm\0 



+ 3HSp,n = 



e + 2He = 

5G{p^m + pa) + SGpA + G{5prn + ^Pa) + G5pA = 
fc2 [G5pA + Pk5G] + a^pr,^Ge = . 



(29) 

(30) 

(31) 
(32) 
(33) 
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Since Eq. (|3ip will be important for the subsequent considerations, let us note that it follows from 
Eq. (j27p upon using the matter conservation equation ([6]). 

We are now ready to derive an important result that applies to all cosmological models of the 
type iv) in section [21 i.e. models with variable A and G in which matter is covariantly conserved, to 
wit: for all these models, the perturbation equations that we have just derived do not depend in fact 
on the wavenumber k. To prove this remarkable result is very easy at this stage of the discussion, 
as it ensues immediately from equations (|3ip and (jSSp . Indeed, if we change the variable from 
cosmic time t to the scale factor a (which is readily done by using f = df /dt = aHdf /da = aHf'), 
equation (f3T]) can be cast as 

+ -0 = . (34) 
a 

Therefore, 

6 = 00 a^^ . (35) 

It follows that the second term on the l.h.s. of Eq. (i33]l behaves as a?pmG9 = 9q PmG. Notice 
that ^0 is the perturbation of the matter velocity at present, which is of course much smaller than 
any value 9i that this variable can take early on, more specifically at any time after the transfer 
function regime has finished (see below). Obviously, 9i = 9Qa~'^ <ti 9q (where Oj <C ao = 1). Thus, 
taking into account that the matter perturbations 6pm are growing rather than decaying, 9 will 
be comparatively negligible. Let us also note that the 0-term in Eq. (|33p is multiplied by the time 
derivative of G, which in all reasonable models (in particular, the one we explore in section [5]) 
should be small. Finally, if we divide Eq. ()33p by k"^ and take into account that in practice we 
are interested in deep sub-horizon scales, i.e. scales A that satisfy Aa ^ (or, equivalently, 
k ^ Hcl)., it follows that the the ^-term in Eq. (j33p is entirely negligible. Thus, in all practical 
respects, this is tantamount to set 9[a) = (Va). The outcome is that the full set of perturbation 
equations becomes independent of fc, as announced^. To better assess the physical significance 
of this important result, we have repeated the calculation in another gauge (the Newtonian or 
longitudinal gauge) and we have obtained the same result for scales well below the horizon (see 
the Appendix at the very end for a summarized presentation and discussion). 

Having shown that the shape of the spectrum is not distorted at late times in our model, 
we may ask ourselves if, in contrast, there are additional sources of wavenumber dependence at 
early times. Recall that the transfer function T[k) parameterizes the important dynamical effects 
that the various fc-modes undergo at the early epochs [35]. More specifically, it accounts for 
the non-trivial evolution of the primordial perturbations through the epochs of horizon crossing 
and radiation/matter transition (the latter occurs at the so-called "equality time", teq)- A most 
important feature in the ACDM case is that the scale dependence is fully encoded in T{k), and 
so the spectrum of the standard model evolves without any further distortion during most of the 
matter epoch till the present time. 

We may ask ourselves if the evolution of our model with running A and G at the early epochs 
follows also the same pattern as the ACDM, or if its dynamics could imprint some significant 
modification on the structure of T{k). In such distortion of the power spectrum with 

^This result is similar to the model of Ref. [21] . where matter is also covariantly conserved. 
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respect to the ACDM would be generated at early times and it would freely propagate (i.e. without 
further modifications) till the present days and act as a kind of signature of the model. However, we 
deem that this possibility is not likely to be the case, for in the kind of models under consideration 
there is no production/decay of matter or radiation. Thus, the value of the scale factor at equality, 
tteq = a{teq), must be identical to that of the standard ACDM model (cf. more details in section 
15. ip . On the other hand, neither the value of pA nor of the perturbation in this variable are 
expected to play an important role in the past. Finally, for reasons that will become apparent 
later, it is important to restrict our discussion to models wherein the time variation of G is very 
small, as this may also affect the previous consideration on the transfer function (and on top of 
this there are important bounds from primordial nucleosynthesis to be satisfied, see section [5^2]) . 

After insuring that these conditions are fulfilled by the admitted class of variable G models, 
we trust that the ^-dependence of the transfer function for these models should be the same as 
that for the ACDM model. Combining this fact with the above proven scale independence of the 
late time evolution of the coupled set of perturbations for 6p\ and 5G, we reasonably infer that 
the matter power spectrum of a model with variable p\ and slowly variable G - and with self- 
conserved matter components ~ must generally have the very same spectral shape as that of the 
ACDM model. Fortunately, in spite of their shape invariance, such models can still be constrained 
and distinguished from the standard cosmological model by means of the spectrum amplitude (i.e. 
from the normalization of the matter fluctuation power spectrum) . We shall further elaborate on 
these points in section [5l where a concrete model exhibiting these properties will be analytically 
and numerically analyzed. 

Let us now retake our analysis of cosmic perturbations by showing that the perturbations in 
G are actually tightly linked to the perturbations in p\. This feature is another consequence of 
the aforementioned scale invariance of the late time evolution of the perturbations. Indeed, from 
([55]) and the neglect of the velocity perturbations ([55]) . we obtain: 

PA G 

It follows that the two kind of perturbations are not ultimately independent. Furthermore, using 
Eq. ()36p . a straightforward calculation from ()32p renders the simple result: 

_6p^ 6G {5Gy 

The last two equations show clearly that the perturbations in G and p\ may not be consistently 
set to zero, since they will be generated even if they are assumed to vanish at some initial instant 
of timelll From (pOl) . and using the matter conservation ([6]), it is easy to see that 

/i = 2 . (38) 

We can now use this last expression and Eq. (|29p to produce a second order differential equation 



*This is analogous to what happens in models with self-conserved DE, where DE perturbations cannot be con- 
sistently neglected [36] 
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for 5m- Using again differentiation with respect to the scale factor, we find 

C + A{a)C = B{a){8m + ^-§^ , (39) 



where we have defined: 



^ = - + ^ (40) 



a H{a) 

3 
2^' 



B = ^hm{a) (41) 



^m{a) = — -— = . Pm{a) . (42) 

Pc[a) 3H^[a) 

Here, Qmict) is the "instantaneous" normahzed matter density at the time where the scale factor 
is a = a{t). It is also convenient to define the corresponding instantaneous normalized CC density 
0,\{a) = PA(a)/pc(o)> and it is easy to check that the following sum rule 

Qmiz) + QAiz) = 1 (43) 

is satisfied for all z. Clearly, the sum rule is a simple but convenient rephrasing of Eq. ([3]). 

If we would neglect the perturbations in G {6G ~ 0) and hence 6pA ~ too, owing to Eq. ([36]) . 
it is immediate to see that Eq. (j39p would reduce to the standard differential equation [35] for the 
growth factor for DE models with self-conserved matter under the assumption of negligible DE 
perturbations, i.e. explicitly Eq.(90) of Ref. |21j . In particular, if we take for H the standard 
form ([3]) with pA =const., we arrive to the equation for the ACDM model. However, Eq. (I39p with 
6G 7^ tells us precisely how to take into account the effect of non-vanishing perturbations in the 
gravitational constant G. We could now solve the coupled system formed by ()39p and (j37p . giving 
initial conditions for 5'^, 5m and 5G, or use those two equations to get a third order differential 
equation for 5m, in which case we would need to give the initial condition for 5m instead. This seems 
to be the natural thing to do, because it involves boundary conditions on 5m and its derivatives 
only, and does not mix with the boundary conditions on other variables. 

In order to get the aforesaid third order differential equation for 5m, we have to differentiate 
Eq. ([39l) and use this same equation to eliminate 5G/G. Using also (f37|l . we obtain 

5G\ 5G' G'5G G' G' f 5l + A5'm .\ G' . 
^ =^-G^ = -G^--g[^ '^^ =-^('^™ + ^'^J, (44) 



and hence we finally arrive at the desired third order equation: 

A' -B + A 



6'Z+iA-^ + ^]5^ + 



G[_B^ 



5'm = 0. (45) 



Now we have to compute the coefficients of this equation. From the definition of Qm, Eq. (|42p . we 
get the useful relation: 

nL G' ^H'_3 ^^g^ 



Qm G Ha 
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On the other hand, differentiating Priedmann's equation ([3]), 

2HH' = ^[G'(p„ + pa) + G{p'^ + p'j,) = ^p'^ = -^pm , (47) 

where in the last two steps we have taken into account the energy-conservation law ([8|) and ([9|) for 
o^m = 3 (matter epoch). Combining this last equation with the definition of flm, (|42p . we get: 

Finally, using ()46p and ()48p we can rewrite our third order differential equation (|45p in terms of 

C + l (l6 - 90„) ^-k + ^(^8-im^ + 3nl-an'Jl^-^ = 0. (49) 

Let us stress that this equation is valid for any model with variable G and pjy in which matter 
is covariantly conserved, i.e. models of the type iv) in section [2l However, a very particular case 
where it should give the correct results is for the ACDM and CDM models because matter is 
conserved and the equation ([8]) is trivially satisfied. For the CDM model, for instance, = 1 
and the equation reduces to: 

^ (Jj 



which has the general solution: 



5m = Cia + C2 + C^a-^'^ . (51) 



Barring the constant and decaying modes, which we can neglect, the relevant solution is the growing 
mode 5m oc a, which is linear in the scale factor, as expected. Note that the setting VLm = 1 implies 
that there is no CC term, and from Eq. ([8|) we infer that G must be strictly constant in such case. 
Finally, Eq. (|36p entails that 5G is then also constant, and this constant can be set to zero through 
a redefinition of G. Let us also remark that, in this limit, one does not expect to recover the 
perturbative analysis of scenario ii) in section [21 i.e. the results obtained in [28]. The reason 
is that although in the last reference there are no perturbations on G, matter and vacuum are 
exchanging energy. Therefore the underlying physics is completely different and there is in general 
no simple connection between these two kind of models. 

We can use the general equation ()49p to study the perturbations in any model with variable G 
and PA in which matter is conserved. In the next section, we consider a well motivated non-trivial 
example. 



4 The class of vacuum power series models in H 

As we have mentioned in the introduction, there have been many attempts to envisage a dynamical 
cosmological term p\{t) [23]. However, in most cases the approach is purely phenomenological, 
with no reference to a potential connection with fundamental physics, via QFT or string theory. 
On the other hand, there are models wherein there is such a possible connection. Obviously, this 
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represents an additional motivation for their study. Consider the class of models in which the CC 
evolves as a power series in the Hubble rate: 



where rii [i = 0,1,2,...) are dimensional coefficients (except 71,4, which is dimensionless) . The 
background solution for this cosmological model up to z = 2 can be found in [TS]. The higher 
order terms in (j52p . made out of powers of H larger than 2, are phenomenologically irrelevant at 
the present time and will not be discussed. Besides, let us stress that from the fundamental point 
of view of QFT in curved space-time, the general covariance of the effective action can only admit 
even powers of the expansion rate |15 y i9 t [57]. So the coefficient ni of the linear term in H cannot 
be related to the properties of the effective action, but to some phenomenological parameter of the 
theory (e.g. the viscosity of the DE fluid) which is not part of the fundamental principles. If we 
dispense with this kind of phenomenological coefficients, we obtain the subclass of models of the 
form 



which have been advocated as scenarios in which the CC evolution can be linked to the renormal- 
ization group (RG) running in QFT |15tll6j. Notice that the coefficient n2 has dimension of an 
effective mass squared, ri2 = M'^^. We shall further comment on it below. 

The form (j53p is indeed crucially different from just considering that the vacuum energy is 
proportional to , in the sense that Eq. (j53|) is an "affine quadratic law" (i.e. no 7^ 0). While 
the pure law is not favored by the experimental data [TH], the affine version has been recently 
tested in a framework where G is constant, specifically in the framework of the scenario ii) in 
section [21 and it was found that it can provide a fit to the current observational data of similar 
quality as the ACDM one - see [18] for details. 

Consider now the size of the 71-2 coefficient in Eq. (j53p . In its absence, the the CC is strictly 
constant and uq just coincides with the current value p^. However, in the presence of the 
correction, the boundary condition at H = Hq becomes = uq + 77-2 -f^o- additional term 

is to play a significant role it should not be negligible as compared to uq, and as the latter is the 
leading term it must still be of order ~ /O^. It follows that the effective mass M^g associated to the 
coefficient n2 cannot be small, but actually very large - specifically, of order of a Grand Unified 
Theory scale (see below) . This also explains why no other even power of H can play any significant 
role in the series expansion (j52p at any stage of the cosmological history below a typical GUT scale. 
The upshot is that, in practice, the evolution law (|53p is the leading form throughout all relevant 
cosmic epochs (radiation dominated, matter dominated and late CC dominated epochs). 

Likewise, the possibility that the vacuum energy could be evolving linearly with H ~ i.e. as if 
riQ = n2 = in Eq. (j52p - has also been addressed in the literature and can be motivated through 
a possible connection of cosmology with the QCD scale of strong interactions [38fl4U| . However, 
as we have said, this option is not what one would expect from the the general covariance of the 
effective action. Actually, the confrontation of the purely linear model H with the data does not 
seem to support it |181l41j. and therefore the linear term alone is unfavored. However, it could 
perhaps enter as a phenomenological term in a general power series vacuum model of the form 



Pa{H) =no + niH + n2H'^ + ... 



(52) 



Pk{H) = no + n2i?^ 



(53) 
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- a possibility which is currently under study [12] • For some alternative recent models with 
variable cosmological parameters, see e.g. |43j . 

In the following, we focus on the quantum field vacuum model of the form (jSSp . It is particularly 
convenient to rewrite the coefficients of this equation as follows: uq = p\ — Mp Hq / (Svr) and 
n2 = Su Mp/{87r). Therefore, the CC evolution law reads 

PaW=pU^^4(^'-^o)- (54) 

Here i/ is a small coefficient (|z/| <^ 1) and Mp is the Planck mass; it defines the current value of 
Newton's constant: Gq = 1/Mp. Clearly, the vacuum energy density (j54p is normalized to the 
present value, i.e. 

PAiHo) =pl^^ni Hi Ml , (55) 

where experimentally ~ 0.7. The above parametrization satisfies the aforementioned condition 
that no is the leading term and is of order p^, and at the same time the correction term is of order 
Mggiif^, with Meff ~ ^/v Mp a large mass, even if v is as small as, say, ~ 10~^ or less. 

It is now convenient to define a new set of cosmological energy densities normalized with respect 
to the current critical density p^^ = 3Hq/{87t Gq): 

n^(z)^Pli^ = ^n,iz) (z = m,A), (56) 

where we have introduced the ratios 

E{z) = ^ = [^m{z) + n^izf^ , g{z) = ^ , (57) 

with G{z) Newton's constant at redshift z. Notice that, in Eq. ()56p . we have explicitly related the 
new parameters ^i{z) with the old ones Q.i{z), the latter being defined in Eqs. (I42p .( l43p . Obviously, 
they all coincide at z = with the normalized current densities il^, i.e. r2j(0) = Oj(0) = VL^. 

Clearly, the coefficient v in Eq. (|54p measures the amount of running of the CC. For any given 
we can compare the value of the dynamical CC term at a cosmic epoch characterized by the 
expansion rate - or by the redshift z - with the current value (|55p . The relative correction can 
be conveniently expressed as follows: 

AnA(z) = M^l^ = [E\z) - 1] , (58) 

Of course G(0) = Gq. Moreover, since g = 1 for = 0, it follows that for small i/, g{z) deviates 
little from 1, namely g{z) = 1 + 0{v). Thus, expanding to order v in the matter epoch, it is easy 
to show from the previous equations that 

A17a(^)=^z.^ [(1 + ^)'-1] , (59) 

where g{z) ~ 1 to this order. If we look back to relatively recent past epochs, e.g. exploring 
redshifts z = 0(1) relevant for Type la supernovae measurements, we see that AQa(-z) can be 
of order of a few times v. For example, for z = 1.5 and z = 2, we have AQa(1-5) ~ Qu and 
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AQ/^(2) ~ 11 respectively, assuming = 0.3. The correction is thus guaranteed to be small, 
as desired, but is not necessarily negligible. We shall see, in the next sections, the potential 
implications on important observables, which will actually put tight bounds on the value of v. 

Although the above parametrization of the CC running can be purely phenomenological, let us 
recall that the dimensionless coefficient v can be interpreted, in more fundamental terms, within 
the context of QFT in curved-space time; specifically it is proportional to the "/3-function" for the 
RG running of the CC term. The predicted value in this QFT framework is |15 y i9 1 [M l l25j 

ly = ^ , (60) 

127r M2 ' ^ ^ 

where M is an effective mass parameter, representing the average mass of the heavy particles of the 
Grand Unified Theory responsible for the CC running through quantum effects (after taking into 
account the multiplicities of the various species of particles). Obviously, M ~ M^s- Since a = ±1 
(depending on whether bosons or fermions dominate in the loop contributions), the coefficient v 
can be positive or negative, but |z^| is naturally predicted to be smaller than one. For instance, 
if GUT fields with masses Mj near Mp do contribute, then < l/(127r) ~ 2.6 x 10"^, but 
we expect it to be even smaller in practice because the usual GUT scales are not that close to 
Mp. By counting particle multiplicities in a typical GUT, a natural estimate lies in the range 
u = 10^5 - 10"3 (see [M1E5] for details). 

In the next section, we shall concentrate on a specific running QFT vacuum model of the type 
(j54p where the vacuum energy and the gravitational constant vary simultaneously in accordance 
to the Bianchi identity We shall also perform a detailed numerical analysis of our results. 



5 Application: running QFT vacuum model with variable G 

In this section, we consider a CC running vacuum model in which p\ = pk{H) depends on the 
Hubble rate through the affine quadratic law (|54p and in which matter is separately conserved - 
the time variation of the CC being compensated by that of the Newton's coupling, G = G{H). 
This setup corresponds to the scenario iv) , as defined in section [2] and was first considered at the 
background level with alternative motivations in Refs. [241125] . The novelty here is to consider 
the detailed treatment of the cosmic perturbations in that scenario. Namely, we apply to it the 
general treatment of cosmic perturbations with variable /oa and G developed in section [3l In this 
way, we can study the growth of matter density perturbations and obtain a LSS bound on the 
basic parameter u, the one that controls the variation of G and p\. Actually, the tightly bounded 
region will appear in the (i^, 0^) plane. 

Specifically, we will constrain the parameter v (j60p by requiring that the amount of growth of 
matter perturbations in our model does not deviate too much from the value deduced from the 
observations of the number density of local clusters. In fact, let us recall that these observations 
allow to set the normalization of the matter power spectrum [44^45) . and hence fix its amplitude. 
In view of the shape independence of the power spectrum for the models under study, which we 
have discussed in section [3l the truly relevant parameter to constraint our model is the spectral 
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amplitude [33]. As we shall see, the constraints obtained in this way are in very good agreement 
with those arising from primordial nucleosynthesis. Throughout this section, we will be using the 
scale factor and the cosmological redshift z = {\ — a)/a interchangeably. 

The cosmological evolution of our model is determined by the Friedmann equation ([3]), the 
Bianchi identity ([SD and the RG law for the CC (|54p . Using the normalized density parameters 
defined in Eq. ()56p . the basic cosmological equations can be formulated as 

E^{z)=g{z)[^m{z) + ^K{z)] , (61) 

{^n, + ^A)dg + gdnK = Q, (62) 

VLf,{z) = ^l + u[E'^{z)-l\ , (63) 

n^{z) = Ql{l + zf^^+^-\ (64) 

where the last equation just reflects the covariant conservation of matter, i.e. it is a rephrasing of 
Eq.Q. Solving the remaining system for the function g = g{H), it is easy to arrive at 

We confirm that g{H) depends also on the parameter and that, to first order, we have g = 
1 + 0{u). Thus, in this model, u plays also the role of the /3-function for the RG running of G. 
Compared to the quadratic running of the CC with the Hubble rate, indicated in Eq. (I63p . the 
running of G with H is indeed very slow, it is just logarithmic. The functions g{z) and ^a{z) 
cannot be determined explicitly in an analytic form. Nevertheless, it is possible to derive an 
implicit equation for g{z) by combining (|65p with (j6ip and ()63p . For the matter-dominated epoch, 
the final result reads: 

1 



^ 1 + In 



9{z) 



9{z) 



In [nm{z)+^l-v\ , (66) 



with nm{z) = n'^{l + zy. The CC density follows from ([63]), and is given as a function of g{z): 

1-u g{z) 

As a simple check, we can see that Eq. (I66p is satisfied at z = for any value of using g{Q) = \ 
and the sum rule + = 1 for flat space. Then from ()67p we immediately obtain r2A(0) = 
as expected. 

In Fig. [H we show the evolution of different background quantities in this model in terms of 
the redshift z. For the plots, we have used = 0.3 and both a positive {v = +4 x 10^^, red 
solid line) and a negative value {v = —4 x 10~^, blue dashed line) for the parameter v. As we will 
see later in this section, according to LSS and nucleosynthesis considerations, we expect u to be 
(at most) of order 10~^, so these are reasonable values. In Fig. [1^, we plot the evolution of the 
matter density fraction, Om(-z) (|42p . which rapidly approaches unity in the past. This means that 
^\{z) is negligible for high redshifts - recall the sum rule (|43p - so that in this regime our model 
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Figure 1: Evolution of different background quantities for the QFT cosmological model with 
running A and G, in terms of the redshift z. We take = 0.3 and consider both a positive 
{y = +4 X 10"^, red solid lines) and a negative value {v = —4 x 10~^, blue dashed lines) for the 
parameter z^, which controls the variation of G and A; a) the matter density fraction, Vlra{z) (j42p . 
rapidly approaches unity in the past, meaning that for high redshifts ^\{z) is negligible and our 
model closely resembles a CDM (for which Om(z) is exactly 1); b) evolution of G{z) (thin lines) 
and H'^{z) (thick lines), normalized to the ACDM values, showing that the departure from the 
standard cosmological model is small. The evolution of both quantities is related by Eq. (|7ip . 

closely resembles a CDM model (for which we would have ^rn{z) exactly 1.) This could have been 
anticipated from (I63p . which far in the past reads: 

nx{z) ^vE^{z) . (68) 

Using this expression in (|56p . we obtain the corresponding asymptotic value of Vt\ in the past: 

nt,{z) = ^^^n^{z)^vg{z)^0{u)^l (z»l), (69) 

since indeed g{z) ~ C(l), as can be confirmed from Fig. [T}3. Equation (|69p tells us that, for 
any reasonable value of i^, the contribution of the running CC in the past will be unimportant. 
Nevertheless, as seen in the detail frame of Fig. [TJi, Vtm{z) in our model is indeed not exactly one, 
nor really constant. Note that for < the CC density decreases as we go into the past (the 
opposite is true for u > until it eventually gets negative. Therefore, our model can accommodate 
either originally positive or negative values for A, thanks to the running nature of this quantity. 
Notice that the previous results can be viewed as being a consequence of the fact that, for flat 
space, the tilded normalized densities satisfy the sum rule (jl3|). 

In Fig. [Ij) we show the evolution of g{z) and H'^{z)/H\{z) where is the (flat) ACDM 
Hubble function, simply given by: 

Hi = Hi [^miz) + ^l] = Hi [Q™(1 + zf + ^ll] . (70) 
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Taking into account (|6ip and (j67|) we obtain, to order u, 

Therefore the evolutions of H'^{z)/H\ and g{z) are expected to be very similar, as indeed shown in 
Fig. [1)3. Furthermore, both quantities stay close to 1, so the deviation from the standard ACDM 
evolution is reasonably small, although it maybe large enough so as to be detected in a future 
generation of precision cosmology experiments. For instance, for the values of u and that we 
are considering, the relative deviations of g{z) and H'^{z) at z = 2 with respect to the 

standard model values are (approximately) 4%, 1% and 0.5%, respectively. 

5.1 Constraints from large-scale structure 

Let us now move to the detailed study of the growth of matter density perturbations in our model, 
and further elaborate on some important issues raised in section [3l The matter power spectrum 
for a sufficiently recent value of the scale factor (a ^ Ocq) can be written as: 

P{k,a) = \6m{k,a)\'^ = Ak''T^{k)D^{k,a) . (72) 

Here ^ is a normalization coefficient; n is the scalar spectral index (which gives us the shape of 
the primordial spectrum; e.g. n = 1 if we assume a Harrison-Zel'dovich spectrum); T{k) is the 
transfer function, which does not depend on the initial conditions but only on the physics of the 
microscopic constituents; i.e. it encodes the modifications of the primordial spectrum that arise 
when taking into account the dynamical properties of the particles and their interactions. The 
T{k) function receives non-trivial contributions only from the evolution of the different (comoving) 
wavelenghts A ~ in the epochs of radiation, horizon crossing and radiation/matter equality, 
i.e. it reffects the /c-dependent features that occur at early times when a moves from a <^ Oeq to 
a ^ Oeq. In general, the form of T{k) depends on the cosmological model under consideration, and 
determines to a large extent the final shape of the spectrum. However, for the late time evolution 
(a ^ CLcq), there might be more distortions of the spectrum in models that depart significantly from 
the standard ACDM scenario (for which there are no further fe-dependence beyond that already 
encoded in the transfer function). These additional, model-dependent, effects are refiected in the 
fe-dependence of the growth factor D{k, a) in Eq. (|72|) . A dependence of this sort would appear e.g. 
in models of type ii) in section [2] where the CC decays into matter, as it was previously studied 
in [28] . However, our claim (formulated in section [3|) is that this is not the case for the variable 
CC and C models under consideration, i.e. models of type iv) in section [2] with self-conservation 
of matter, provided that G is a slowly varying function of time. Obviously this is so for the model 
studied in the previous section, where G varies logarithmically with the expansion of the universe, 
see Eq. (j65p . For these models, and of course also for the standard ACDM model, the growth factor 
is independent of k and can be written as 

Dia) = , (73) 

where 5cDM(flo) is the present matter density contrast in a pure cold dark matter (CDM) scenario, 
taken as a fiducial model. 
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5.1.1 More on the transfer function for variable p\ and G models. 

A few additional comments on the transfer function are now in order. As previously commented, 
we expect the transfer function in our model to be the same as for the ACDM. First of all, let us 
try to better justify this expectation, which implies that the shape of the matter power spectrum 
in our model coincides with that of a ACDM with the same value of fi^. Then we show that 
we can still constrain our model by means of the growth factor, obtaining bounds in very good 
agreement with those arising from primordial nucleosynthesis. 

The scale dependence of the transfer function for CDM models is different for large scales (small 
fc's), where T{k) = 1, and small scales, where it asymptotes to ln(A;)/A:^ [35]. The turnaround 
occurs at the value k = k^q, corresponding exactly to the scale that enters the Hubble horizon 
{H~^) at the moment of matter-radiation equality, i.e. at the point a = Oeq which satisfies 

^miaeq) = ^^r(«eg), henCC 

«eq = 7^ , (74) 

where ~ 10~^ is the present density of radiation. The modes that enter the horizon at 
the equality time, with comoving wavelength Aeq, have a physical wavelength that follows upon 
multiplication with the scale factor ([7¥|) . i.e. Aeq aeq = 1/H{acq), or, equivalently, 

keq = aeqH{aeq) ■ (75) 

Obviously, since H decreases with time the perturbations with wavelenght shorter than Agq (i.e. 
those with with k > A;cq) will enter the Hubble horizon before the matter-radiation equality, i.e. in 
the radiation era {t < teq)- From this moment until t = teq, the growth of inhomogeneities in the 
cold DM component becomes suppressed because the expansion rate during the radiation epoch is 
faster than the characteristic collapsing rate of the CDM. As a result, the modes in the radiation 
epoch can grow at most logarithmically with the scale factor. Only after the cold component begins 
to dominate {t > teq) the amplitude of the formerly inhibited modes starts growing linearly with 
the scale factor. Finally, as light can only cross regions smaller than the horizon, the suppression in 
the radiation epoch does not affect the large-scale perturbations {k ^ k^q), which enter the horizon 
in the matter epoch. Such different behavior of the perturbations according to their entrance in 
the horizon before or after the time of equality is the origin of the characteristic shape of the 
transfer function for CDM-like models in the various available parameterizations [35] . 

From the previous standard discussion, it is apparent that the shape of the transfer function 
depends critically on the value of fccq, which in turn depends on Oeq and H{aeq)- In the models we 
are studying, dark matter and radiation are separately conserved, and therefore the value of a^q 
will not change with respect to the standard one, Eq. ()74p . So the remaining issue is to clarify if 
the change of H{aeq) in our model as compared to the corresponding ACDM value is significant 
or not. The answer follows easily from Eq. (j71|) . At the high redshift where equality of matter and 
radiation occurs, z = 0(10^) » 1, the function that accompanies u on the r.h.s. of that equation 
is virtually equal to one, and we are left with 



F(aeq) ~ \ giaeq) (1 + 1^) /?A(aeq) . (76) 



18 



Here, as before, H\ is just the standard ACDM Hubble rate. For the maximal values of u that we 
will be considering (|z^| ~ O(10~^), see the subsequent sections), it is easy to see from ([65]) that 

,|^| ln:%al^l%, (77) 



g{a. 



cq J 



where the expansion rate at the time of equality is //(a^q) ~ 10^ Hq. We see that the change of 
Hi 

d-eq) with respect to H\[aQq) is mainly caused by the variation of g from 1 at t — teq- However, 
numerically, the effect is very small. These results allow us to safely conclude that the value of kcq 
expected in our model is essentially the same as that of the ACDM model, up to differences of 1% 
at most. Given that in the ACDM model the (comoving) wavenumber at equality is 



feeq = acqi?A(acq) = W ^ ^m^O ~ 10"^ h Mpc'^ , (78) 







we see from this expression that a 1% change in 5 (acq) ~ hence in H{acq) - is equivalent to a 1% 
change in Q,^^ in the standard scenario (i.e. with g = 1). This variation is too small to be within 
reach of the present observations (see the next section), and can be safely neglected. 

The final point is that, as argued throughout this section, neither p\ nor its perturbations or 
those of G are important in the past. Therefore, the evolution of the perturbations both in the 
radiation and in the matter-dominated eras (and both in the case of sub-Hubble and super-Hubble 
perturbations) remains essentially unchanged. 

All in all, we expect that the transfer function in the model under consideration is very close 
to that of a ACDM with the same value of From the line of our argumentation that we have 
used, it is not difficulty to convince oneself that this result can be extended to any model with 
variable A and G, and self-conserved matter, as long as G{a) is not changing too fast. Adding 
this property to the fact - cf . section [3] - that the late growth of matter density perturbations in 
these models does not depend on the wavenumber, the final robust conclusion is that the matter 
power spectrum for models within the scenario iv) of section [2] will present the same shape as 
in the standard ACDM case. Therefore, the spectrum shape will not be useful to constrain the 
additional free parameters of the model. This is in sharp contrast to models within scenario ii), 
in which there is an exchange between dark matter and vacuum energy. For these models there 
is an explicit dependence on k beyond that of the transfer function and this produces a late time 
distortion of the power spectrum with respect to the ACDM. As indicated, this was exemplified 
in Ref. |28y29| for the case of an evolution law of the type (|53p . 

To summarize, there is a significant difference between the running QFT vacuum model (|53p 
when studied either in scenario ii) or when considered in scenario iv). While in Ref. |28] I29| the 
shape of the spectrum was used to restrict the parameter v for type ii) models, in the next section 
we show that for the alternative type iv) models, being them shape-invariant with respect to the 
ACDM model, one can make use of the amplitude of the power spectrum in order to constrain the 
free parameters. 
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5.1.2 The spectrum amplitude as a way to constrain running G and A models. 

The normahzation of the matter power spectrum on scales relevant to large-scale structure can 
be performed through different methods [IB], e.g. from the microwave-background anisotropics or 
by measuring the local variance of galaxy counts within certain volumes. One of the most robust 
ways to do it is through the number density of rich galaxy clusters |44y 45| . which is very sensitive 
to the amplitude of the dark matter fluctuations that collapsed to form them. The typical scales 
for these fluctuations are of order 10 Mpc, which are the smallest ones still in the linear part of 
the spectrum. The cluster method determines the amplitude of the power spectrum on just that 
length scale (corresponding to wavenumbers in the upper end of those explored with this method). 
The normalization is usually phrased in terms of as, the root mean square mass fluctuation in 
spheres with radius 8h~^ Mpc (i.e. ~ 10 Mpc, for h ~ 0.7). By assuming a ACDM model, these 
studies are able to provide constraints in the ag - plane. An important feature is that the 
results are approximately independent of the spectrum shape and galaxy bias [H]. For instance, 
in this last reference, it is found the relation: 

as = (0.495l°:°rr) {^'J-'-'' , (79) 

valid for spatially flat models with a wide range of shapes (in particular, valid for 0.2 < Q,'^ < 0.8). 
When combined with CMB analyses, cluster studies can determine both cjg and fi^. In a recent 
work [35], local cluster counts are used in conjunction with WMAP5 data to find: 

0^ = 0.30;H3 (68%). (80) 

We will use this result to constrain our G-variable model. In order to do this, we first define the 
"amount of growth" , namely the square of the growth factor ([73]) at present, /^^(ao), which appears 
directly in the formula of the power spectrum, Eq. ()72p . The idea is to compare the amount of 
growth in our model model with the amount of growth in the ACDM model with Q,^ = 0.30. 
From the standard expression for the growth factor in the ACDM model [35] one finds that a 
~ 10% variation in given by (|80p represents a ~ 5% change in the amount of growth. As the 
determination ([80]) entails only a la constraint, we will be conservative and allow up to a 10% 
deviation in the amount of growth of our model with respect to the ACDM model. We are thus 
asking our model to pass the following "F-test" [291137] : 



F 



i^2(ao, 0.3,0) 



<0.1. (81) 



Notice that -D^(ao, 0.3, 0) in the denominator is just the amount of growth in the ACDM for the 
central value of (j80p , whereas in the numerator we have the amount of growth for the model under 
consideration at a given non- vanishing value of the relevant parameter u, with left as a free 
parameter. As a result, the constraint (j8ip will generate contours in the u - plane which will 
define the allowed region in parameter space. 

The admissible values for ilJri should, in principle, be those compatible with the shape of the 
spectrum. However, the shape has been measured by the 2dFGRS and SDSS surveys, existing 
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a significant difference between their results. On the one hand (assuming a Hubble parameter 
h = 0.72), the SDSS main galaxy analysis @] favored the result 

= 0.296 ± 0.032 . (82) 

Similar values around ~ 0.3 were found by alternative analyses of the SDSS catalogue [i8] . 
On the other hand, the 2dFGRS collaboration [3] found a much lower matter density 

= 0.231 ± 0.021 , (83) 

obtained from measurements of clustering of blue-selected galaxies. The inclusion of Luminous 
Red Galaxies (LRGs) in the SDSS analysis seems to even increase the discrepancy. For instance, 
the authors of [l9] find = 0.32 it 0.01, although a lower density is recovered when restricting the 
analysis to large scales (1)^ = 0.22 ±0.04 for 0.01 < A: < 0.06 /iMpc"^) It is widely believed that 
the differences are due to the scale-dependence of the galaxy bias (which is apparently stronger 
for the red galaxies that dominate the SDSS catalogue) or even to some kind of systematic effect, 
but the problem has not been fully settled so far [50]. As an example, the last result by the 
SDSS team [5], obtained from the analysis of a LRGs sample (in combination with WMAP5 data), 
yields = 0.289 ±0.019, stiU much larger than (and incompatible with) the 2dfGRS result, ([83]) . 
although if we make allowance for a 10% gaussian uncertainty in h in the 2dfGRS analysis can 
bring both results within la. 

At the end of the day, and taking into account the results from 2dFGRs and SDSS, we think 
that it would be premature to discard any value for Q,^^ in the range (0.21, 0.33) on the grounds 
of structure formation data - as long as it predicts the right amount of growth, e.g. by satisfying 
the F-test (|8ip . Therefore, in our analysis, for illustrative purposes, we will consider values of Q,^^ 
between 0.2 and 0.4. 

5.1.3 Numerical results. 

In this numerical section, and in view of the previous considerations, we wish to determine the set 
of points in the - Q,^ plane for which the amount of growth, determined by the value of D^(ao), 
in our running p\ and G model deviates less than 10% from the central (O^^ = 0.3) ACDM value 
in Eq. ()80p . To compute the growth factor in our model, we evolve the solution 5m (o) of the 
differential equation (|49p from an initial value a = up to the present moment (qq = 1)5 where 

^ 1 is the scale factor at some early time, deep into the matter-dominated era but well after 
recombination (so that the transfer function regime has already ended). We will take Oj = 1/501 
(i.e. Zi = 500), although we have checked that the results do not the depend significantly on the 
specific value. As we have seen in Fig. [1^, early on at a = our model is very similar to the 
CDM model {Qm — 1)) for which the matter density perturbations grow linearly with the scale 
factor, i.e. D{a) = a. Therefore, we will assume that this is also the case for our running G and 
p\ model, and take D{ai) = Oj, D'{ai) = 1, D"[ai) = as the initial conditions for the third-order 
differential equation (jl9|) . 

The results are shown in Fig. [2^. The shaded areas represent the points allowed by our analysis. 
Specifically, we compare the case with perturbations in p\ and G (green band) with the case in 
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Figure 2: Analysis of the parameter space of the running QFT models. The shaded areas represent 
the points for which the amount of growth (D^(ao)) of the model deviates less than 10% from the 
ACDM value. For the latter, we take = 0.3, Eq. (I80p . on the basis of data on the matter 
power spectrum amplitude; a) Scenario iv) section [2] (running A and G, self-conserved matter). 
The green (resp. yellow) band include (exclude) perturbations in p\ and G. With perturbations, 
the constraint on the parameter space becomes tighter; b) Comparison between two scenarios of 
section [2] when perturbations in p\ and G are neglected: scenario iv) (yellow band) and scenario 
ii) (orange band). The deviations from the ACDM case are larger in scenario ii) owing to the 
production/decay of matter from the running pA at fixed G, and the constraint is correspondingly 
tighter. The blue and red points in the plane have coordinates {v = ±4 x 10~^, $7^ = 0.3) - used 
in Figs. [Hand [31 The dashed horizontal lines signal the la limits on from Eq. ()80p . 

which these perturbations are neglected (yellow band). In the last situation, the growth factor 
is obtained by solving Eq. (f39ll . although, as discussed in section [3l it is not possible to neglect 
them consistently. The most remarkable conclusion that emerges from this numerical analysis is 
that by considering the effect of the perturbations results in narrower domains, meaning that the 
deviations with respect to the standard ACDM amount of growth tend to be larger - which is 
why the restriction is accordingly tighter. Therefore, the outcome of the analysis with 5G 7^ is 
that, for any in the range 0.2 < < 0.4, values of \v\ larger than 10~^ are ruled out by 
the data on the number density of local clusters. In particular, this excludes the canonical value 
\u\ = l/127r ~ 0.026 obtained from the simplest choice M = Mp in Eq. (|60p. Notice that for 
in the narrower (Icr) interval (I80p . the allowed values for \v\ are of order 10~^ at most. For 
these values of z/, our model is on equal footing with the ACDM as far as structure formation is 
concerned. 

In Fig. [2)3 we compare the results for two of the scenarios of section [2j scenario iv) , represented 
by the yellow band (this is the scenario we have been analyzing so far, with variable pA and G and 
self-conserved matter; and scenario ii), indicated by the orange band; for the latter, G is constant 
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and the CC exchanges energy with matter. In both cases, we are neglecting the perturbationqj in 
PA, which in scenario iv) imphes that 5G = as well, cf. Eq. (j36p . The effective equation for the 
matter density contrast in scenario ii) is the following: 



* d f ^ 

^T^Gpm - 2H — 

Pm at \pm 



'•m = 0, (84) 
^ = Pm + 3Hp^ = -PA . (85) 



This equation, whose primary derivation was performed within the Newtonian formalism [51] . 
can also can be derived from the general relativistic treatment of perturbations [21j . as explained 
in [18]. It is convenient to express it in terms of the scale factor a: 



a H pmtia^ 



3 ~ 2 ^ a ( 



% = 0. (86) 



2 H pm H \p, 

For = 0, we recover equation ([39|) (with 6G = 0), as expected. In Fig. [2)3, we see that the 
differences in the amount of growth with respect to the ACDM case are larger for scenario ii); the 
natural interpretation is that this is caused by the production/decay of matter associated to the 
time evolution of p/\_. 

The blue and red points in Figs. [2^,b correspond to the values analyzed in Fig. [H i.e. = 0.3 
and u = ±4 x 10~^. We will use again these values to exemplify the evolution of the matter and 
Newton's coupling perturbations. Fig. [3^^ shows the growth factor as a function of the redshift, 
both when allowing for perturbations in G and p\ and when they are neglected. For < 0, our 
model predicts more growth than the ACDM model. When 5G = 0, this is just due to the fact that 
p\ is decreasing when we go into the past, so its repulsive effect diminishes. When considering the 
perturbations in G and A, the deviation with respect to the standard case increases even more, 
as already commented in Fig. [2^. The opposite situation occurs for > 0; here the suppression 
of growth is attributed to the enhanced repulsion of matter associated to a positive pA, which 
is increasing with z. Such suppression is larger when we include the G and A perturbations. A 
detailed study of these effects for model ii) of section [2] was performed in Ref. [22], in an effective 
approach with no perturbations in the CC. 

In Fig. [3)3, the evolution of 5G/G, as obtained from Eq. (j37|) . is depicted. This equation only 
determines 6G' , so we need to give an initial condition for SG at a = a^. In order to do so, we note 
that the fact that Clm{z) ~ 1 in the past (reflected in Fig. [T^) ensures that the perturbations in G 
are not playing any important role, according to our discussion at the end of section [3l Therefore, 
the most natural condition seems to be 6G{ai) = 0, and, besides, we expect 6G/G to remain 
negligible until very recent times {z < 10), when VLrn begins to depart from 1. This is indeed what 
can be seen in Fig. giving additional support to our assumption. Let us remind from (|36|) that 
5a = —5G/G. Thus, for < we have 6G > and hence 6\ < 0, which explains the further 
enhancement in the growth of matter perturbations with respect to the case with 6G = 5p\ = 0. 

The main conclusion of this section is that for v of order 10~^ or smaller, our model is in 
perfect agreement with recent data on the normalization of the power spectrum. For the time 



^Remember that within scenario ii), including the perturbations in A causes the growth factor to become scale- 
dependent [28], D = D{k, a), so the kind of analysis we are performing here would be no longer possible. 
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Figure 3: Evolution of the perturbations for the QFT model with running A and G, using the 
same values of 0,^ and v as in Fig. [U a) The growth factor D(z), showing an enhancement 
(suppression) of the growth with respect to the ACDM case for negative (positive) i'. The effect is 
present even when we neglect the perturbations in G and p\, but it is larger if we consider them; 
b) Evolution of the perturbations in G, under the natural assumption that they are negligible 
at Zi = 500. The perturbations remain unimportant until very recent times, when begins to 
depart significantly from 1 (cf. Fig. [1]). Remembering that 5\ = —5G/G, we can explain the 
enhancement (suppression) of matter perturbations observed in a) on account of the opposite sign 
of 5a = (5pa/pa, see Eq. ([39]). 



being, there is still some controversy on the value of resulting from the power spectrum shape 
measured by the two main galaxy surveys. For a given value of il^, the power spectrum in our 
model presents, in very good approximation, the same shape as in the ACDM model. However, a 
non-vanishing ly could become manifest through a difference in the amount of growth. Therefore, a 
future simultaneous precision cosmology measurement of the shape and amplitude (normalization) 
of the matter power spectrum should have the potential to discriminate between a model of variable 
PA and G, and the ACDM. 

In the following section, we show that u ~ 

0(10-3) is also the maximum allowed value that 
we would expect from arguments related to primordial nucleosynthesis, fact that reinforces the 
validity of this bound. 

5.2 Constraints from primordial nucleosynthesis 

Big Bang nucleosynthesis (BBN) can also provide limits on the possible variation of Newton's 
coupling G. The BBN predictions for the light element abundances are sensitive to a number of 
parameters, such as the baryon-to-photon ratio t] = ns/n^ (in view of the fact that the nuclear 
reaction rates depend on the nucleon density) or the value of the Hubble function at the nucle- 
osynthesis time, H]^ = H{z = z^) (given that the expansion rate competes against the reaction 
rates). Since rj oc 0,b can be accurately determined through CMB measurements, one can use the 
observed abundances to constrain the expansion rate, and since H oc VG these constraints can be 
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directly translatecifl into bounds on gj^ = g{zj^), where g{z) was defined in Eq. ([^7|) . 

The current constraints on g^ available in the literature usually make use either of the deu- 
terium abundance (D/H) or the ^He mass fraction (Yp). Deuterium has the advantage that it is 
not produced in significant quantities after nucleosynthesis; its primordial abundance can be deter- 
mined in a quite precise manner by studying the spectrum of the light from distant quasars, which 
exhibits an absorption line corresponding to the deuterium present in (high-redshift) intervening 
neutral hydrogen systems. This turns deuterium into an excellent probe of the universe at the 
time of BBN. However, its predicted abundance depends much more strongly on r/ than on Hjsf. 
As a consequence, the constraints on g^ based on deuterium measurements are not very stringent. 
For instance, in [52] it is found g^ = l.OllQ^g at the 68.3% confidence level or: 

gN = lMt'oil (95%). (87) 

The abundance of ^He is much more sensitive to the expansion rate, but accurate observational 
values are difficult to obtain, since there are many potential sources of systematic uncertainties [53]. 
As a result, a wide range of values for Yp can be found in the literature, see e.g. Table 12 in [SI]. In 
this reference, the authors used Yp = 0.250 it 0.003 (in fact, they use Yp in combination with D/H, 
although the dominant effect is that of helium [5^) to obtain (see section 8.2.5 of the latter): 

0.964 <gN < 1-086 (95%) , (88) 

whereas in [55| a more conservative range for the ^He mass fraction is adopted (Yp = 0.2495 it 
0.0092), leading to: 

0.9 < < 1.13 (68%) . (89) 

Applying (j87p - (|89p to our model, we can effectively constrain our parameter v. In order to do so, 
we compute g^ from (j65]) : 

_ 1 1 

where 0^ is the radiation energy density fraction at present, we have neglected both the matter 
and CC contributions to the expansion rate at z = zj^, and is evident that g^ can be neglected as 
well in the logarithm. Taking Q,^ ~ 10~^ (which includes photons and three light neutrino species) 
together with zn ~ 10^, and comparing the last expression to (|87l) - (f89]) . we find: 

. 1.. 1 A 1 1 f 



From dSTD — > -4.1 x 10"^ < < 5.5 x 10 
From dSHD — > -1.1 X 10"^ < < 5.1 x 10 
From (1891) — > -1.6 x 10"^ < < 1.5 x 10 



Let us remind that (I88p was derived from a value for Yp with a (possibly unrealistic) very small 
error and that (I89p is a 68% value (whereas the other two limits are given at the 95% confidence 
level.) In any case, the conclusion arising from this nucleosynthesis analysis seems to be that the 



^Provided that the total energy density at z = zn is approximately the same as in the standard case. Eq. (|69p 
shows that this condition is satisfied in our model. 
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parameter v can be, at most, of order 10" . This is in complete agreement with the constraint on v 
obtained in section ISTT] from structure formation considerations. Therefore, we have arrived to the 
same result by using two very different methods, which gives additional credit to our conclusions. 

6 Conclusions 

In this paper, we have derived the general set of cosmological perturbation equations for FLRW 
models with variable cosmological parameters /9a and G in which matter is covariantly conserved. 
To our knowledge, this is the first time that a complete set of coupled differential equations for bpt\ 
and bG is presented in the literature. We have shown that the linear growth of matter perturbations 
of this model, D{a) oc bm{o)^ is independent of the wave number k. Adding this property to the 
fact that we generally expect these models to present the same transfer function as the ACDM, 
the scaling dependence and hence the shape of the power spectrum will not change with the late 
time evolution and it will coincide with that of the ACDM model. This fact is remarkable as it is 
in contrast to the situation, more frequently studied in the literature, in which the time-evolving 
vacuum models exchange energy with matter at fixed Newton's coupling G |18j . 

We have exemplified the difference in the power spectrum of running vacuum models, with and 
without matter conservation, by considering the class of cosmological models characterized by the 
running CC law (j53p . which we call quantum field vacuum models because the evolution of the 
CC in them is of the form that one would naturally expect from QFT, and more specifically from 
the renormalization group evolution of the cosmological parameters. Such quantum field vacuum 
models depend on a single parameter i^, which acts as the /3-function for the running of the vacuum 
energy density /9a = Pk{H) and the running gravitational coupling G = G{H). It is interesting to 
note that while the running of the vacuum energy pk{H) is quadratical in the expansion rate, the 
running of G{H) is logarithmic in H, Eq. (j65p . and therefore much milder. For v > 0, the running 
of G is asymptotically free, hence G decreases (resp. increases) logarithmically with the redshift 
(resp. with the expansion), whereas for < it decreases with the expansion. The background 
evolution of these models was studied in [MlllS] . 

In the present paper, we have concentrated on the study of the cosmic perturbation equations 
of the running p\{H) and G{H) models. After confronting the predicted matter growth with the 
various cosmological data, we have found that the final region of the parameter space is a natu- 
ralness region in which u could be, in absolute value, of order 10~^ at most. We have also checked 
that this bound is compatible with the restriction on G variation that follows from primordial 
nucleosynthesis. Remarkably enough, the results emerging from the perturbations analysis are 
derived from the amplitude of the power spectrum rather than from its shape. Numerically, the 
obtained upper bound on v is perfectly compatible with the quantum field theoretical definition 
of this parameter, see Eq. ([60]). and it implies that the mass scales that would enter the quantum 
running of the cosmological parameters could be one order of magnitude below the Planck scale, at 
most. The result \y\ < 0(10"^) is also compatible with previous analyses, using various indepen- 
dent procedures, of models with running p\ in which the vacuum can decay into matter jl8 p28f[50] . 
This decay feature, however, is impossible for the models with running pa{H) and G{H) that we 
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have studied here, and this is the basic reason why they can exhibit the same power spectrum 
profile as the ACDM. We cannot exclude that this property could have been responsible for the 
fact that we have observationally missed this fundamental possibility up to now. As we have 
emphasized, we expect that these models should eventually be testable in the next generation of 
precision cosmology observations from the analysis of the spectrum amplitude, rather than of the 
spectral shape. 
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7 Appendix: cosmic perturbations in the Newtonian gauge 

In this appendix, we sketch the derivation of the perturbation equations in the Newtonian or 
longitudinal gauge [56]. Our aim is to explicitly check that the same physical consequences that 
we have derived in section [3] for the synchronous gauge are recovered in another gauge. In this 
way, we confirm in our particular context the general expectation that calculations of cosmic 
perturbations at deep sub-horizon scales should not present significant gauge dependence. 

The most general perturbation of the spatially flat FLRW metric can be conveniently written 
as follows [56] 

ds^ = a^(?7)[(l + 2^^)^ - ujidr]dx' - ((1 - 2(p)6ij + Xij) dx'dx^] , (91) 

where rj is the conformal time, defined through dr] = dt/a. The above metric is expressed in the 
notation of [21] and consists of the 10 degrees of freedom associated to the two scalar functions tp, (p, 
the three components of the vector function ojj (i = 1, 2, 3), and the five components of the traceless 
second-rank tensor Xij ■ Glearly, the synchronous gauge that we used in section [3] is obtained by 
setting ip = 0, oji = and absorbing the function <j) into the trace of Xiji which then contributes 
six degrees of freedom. In this gauge, the metric part of the analysis of cosmic perturbations is 
tracked by the nonvanishing trace of the metric disturbance: h = g^'^ h^y = g^^ hij = —ha /a?. 
The corresponding equations have been presented in detail in section [3] after defining the "hat 
variable": h = -dh/dt = d [ha/ a?) /dt. 

An alternative possibility is to use the (conformal) Newtonian gauge |56y57j. In this case, we 
set bJi = Xij = 0(Vi,j) in Eq. (f9T]) . A useful feature of this gauge is that the metric is diagonal 
and another is that, in the Newtonian limit, ip plays the role of the gravitational potential. It is 
now straightforward to repeat the calculation of the perturbation equations in a similar way as in 
section [3j In the absence of anisotropic shear stress (which is always the case with non-relativistic 
particles, such as baryons and dark matter), the perturbed Einstein's equations for i / j give 
tp = (p [SZ]) see also [58]. Indicating by a prime the derivatives with respect to the scale factor 
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(/' = df /da), the final perturbation equations in that gauge read as follows: 

„2 



k 



5pm + 5pK + SHpmj^O 



= -47ra2|G 



+ [Pm + Pa)SG 



+ 



SG'ipm + pa) + 5Gp'^ + G'i6p^ + 6 pa) + Gdp'^ = 
Ha? ^ 

G6pA + paSG + ' 



(92) 



A;2 



PmG'e = o. 



Notice that the role of h in the synchronous gauge is taken up here by the variable <j) = ip, which 
turns out to satisfy an algebraic rather than a differential equation. Worth noticing is also the fact 
that the last two equations of (j92]) are identical to equations (j32|) and ([33|) respectively, except that 
in the present case we used derivatives with respect to the scale factor rather than with respect to 
the cosmic time (as in this way length scales can be better compared with the horizon). 

We are now ready to show that for deep sub-Hubble (sub-horizon) perturbations, i.e. those 
satisfying k/{aH) ^ 1, these equations lead to the same second and third-order differential equa- 
tions for the normalized matter overdensity 5m = 5pml Pm^ which we have previously found within 
the synchronous gauge in section [3l To start with, the first equation of (j92p can be cast as: 

3 H'^a^ 



2 fc2 



bG_ 
~G 



(93) 



where 5a = 5p\/ p\, and of course we use the same notation as in section [3] concerning the 
parameters (im and (l^. In the previous equation, we have dropped the term 3Hpm a'^ 6/k'^ since it 
is negligible for sub-Hubble perturbations. In this same regime of perturbations, Eq. (l93p obviously 
implies that (j) <^ 5m- Next, let us note that the second equation of (j92p can be simplified by using 
the conservation of matter, giving: 

C + -^ = ^<P'- (94) 
an 

Differentiating this equation and using the third one in (|92p to substitute for 6' , we find: 



5" 



'3 H' 
aH Ka^li 



30" 



k 



,2 ■ 



(95) 



Furthermore, if we eliminate 9 from this equation by means of (j94p and also using (|93p . we get: 

51 - 30" + {5'm - 300 + = ^ [^rnSm + f^A^A + ^} ■ (96) 

Again we can neglect the terms that are small at sub-Hubble scales; for example, in the previous 
equation we can neglect 0' (0") as compared to 5'm {5'/^) on account of the previously discussed 
relation (j) <^ 5m- In this way, we arrive at: 



(97) 
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The almost final step is to substitute 6\ by means of the last equation of (|92p . Once more, we 
can argue (as we did in section [3|) that apart from the damping factor Ha/k <^ 1 for the sub- 
Hubble modes, there are additional suppression effects, such as the presence of G' which is small 
as compared to G, and hence {G' /G)6 can be considered almost second-order. Therefore, the last 
equation of boils down in practice to G6pA + p\5G = 0, which simply renders 5\ = —6G/G. 
Using also the sum rule (143^ . the relation (j97p finally takes on the simpler form 



This equation exactly coincides with Eq. (|39p . showing indeed that our result was not gauge- 
dependent. Obtaining the third-order equation for 6m is now straightforward. From the last two 
equations of ([92]) one can show that 

From here we proceed as in section [3] for the synchronous gauge, i.e. we differentiate Eq. (j98p 
with respect to the scale factor and we reach once more the third order differential equation ()49p . 
(q.e.d.) 

Some discussion is now in order. Let us first note that the variable cosmological constant 
scenario ii) of section [21 whose perturbation equations were originally analyzed in |28j in the 
synchronous gauge, has been reanalyzed in the Newtonian gauge in [32] and the same conclusions 
have been attained. We should also mention a very recent paper [59] where a nice discussion is 
made on comparing the synchronous and Newtonian gauges in the context of a model with self- 
conserved DE and without DE perturbations. As it is well-known, in that case the two gauges 
give the same evolution for the cosmological perturbations for scales well below the horizon - see 
also |57[I60|. This feature was expected in our model too, even if it is a bit more complicated 
than the one considered in (our DE is also self-conserved, but we do have both DE and G 
perturbations), and in fact the physical discussion about the gauge differences is very similar. In 
the same reference, it is also pointed out that the /c-dependent effects are negligible for large values 
of k, although they could start to be appreciable for k < O.Olh (or, equivalently, for length scales 
larger than £ ~ ~ lOOh'^ Mpc). 

The issues about gauge dependence are of course important in general, because scale-dependent 
effects may change from gauge to gauge. For instance, let us note that there is no counterpart 
(within the synchronous gauge analysis considered in section [3]) of the scale-dependent effects 
introduced by the first equation ([^2]) . At the same time, there may be scale dependent features 
which are tied to the particular cosmological model under consideration, and they need not show 
up in another model, even if we compare them within the same gauge. Indeed, take once more the 
model discussed in Ref. [59], where the DE is self-conserved and there are no perturbations in this 
variable. In such situation, the cosmological perturbations of matter in the synchronous gauge turn 
out to be exactly scale-independent (without any approximation). In our case, in contrast, there is 
some residual scale dependence in the synchronous gauge, although it can be shown to be negligible 
- as we have extensively discussed in section [3| Gauge differences, however, can be significant in 
certain cases, specially for calculations involving very large scales, showing that in general it is 
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non-trivial to relate the perturbation variables to observable quantities in different gauges |57y6Uj. 
In that case, one must keep the appropriate ^-dependence (as e.g. in the Newtonian gauge [59]) 
or resort to a gauge-invariant formalism |6HI62j. However, for the range of wavenumbers usually 
explored in the analysis of the matter power spectrum (viz. O.Ol/i Mpc~^ < k < 0.2h Mpc~^ [3]), 
the gauge differences, and in particular the fe-dependence, should be small enough so as to be 
licitly neglected. In other words, for such scales, the sub-Hubble approximation that we have 
taken in this Appendix holds good and the calculations in the two gauges are found to be exactly 
coincident, as we expected. 

Finally, let us note that the reason why the physical discussion can be more transparent in 
the Newtonian gauge is because the time slicing in this gauge respects the isotropic expansion 
of the background. The synchronous gauge, instead, corresponds to free falling observers at all 
points, entailing that its predictions are relevant only to length scales significantly smaller than 
the horizon, but at these scales it renders the same physics as the Newtonian gauge [STJESEni- 
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